BU-HEP/92-17 
TAUP 2034-93 
hep-pri/9304"2T7 
April 1993 



STATIC PROPERTIES OF QUARK SOLITONS 

'&) : 

^ ; Gonen Gomelski* and Marek Karliner^ 

School of Physics and Astronomy 
<^ \ Raymond and Beverly Sackler Faculty of Exact Sciences 



in 



Tel-Aviv University, 69978 Tel-Aviv, Israel 



> \ Stephen B. Selipsky 1 

o 

^ . Department of Physics, Boston University 

O ! Boston, MA 02215-2507, U.S.A. 

m 

On 
43 



Abstract 



It has been conjectured that at distances smaller than the confinement scale but 
large enough to allow for nonperturbative effects, QCD is described by an effective 
SU(N c xN f ) L x SU(N c xN f ) R chiral Lagrangian. The soliton solutions of such a La- 
grangian are extended objects with spin i. For iV c =3, Nf=3 they are triplets of color 
and flavor and have baryon number |, to be identified as constituent quarks. We inves- 
tigate in detail the static properties of such constituent-quark solitons for the simplest 
case Nf=l,N c =3. The mass of these objects comes from the energy of the static 
soliton and from quantum effects, described semiclassically by rotation of collective 
coordinates around the classical solution. The quantum corrections tend to be large, 
but can be controlled by exploring the Lagrangian's parameter space so as to maximize 
the inertia tensor. We comment on the acceptable parameter space and discuss the 
model's further predictive power. 
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1 Introduction 



One of the outstanding questions in particle physics continues to be the relationship between 
the phenomenologically successful non-relativistic constituent-quark model of hadrons, and 
QCD's fundamental degrees of freedom The u, d, s valence quarks appearing in 
hadron spectroscopy have the same quantum numbers as the QCD current-quark fields, yet 
are much heavier and in some sense embody the strongly interacting QCD dynamics. The 
experimental results || on spin structure of the proton indicate that the constituent quarks 
ought to be thought of as composite objects with internal structure. Finding a theoretical 
description of this structure is an important and interesting challenge. 

Kaplan has put forward a model 0] combining some of the features of the chiral quark jl] 
and skyrmion|| ^ [L0| approaches. In this model it is postulated that at distances smaller 
than the confinement scale but large enough to allow for nonperturbative phenomena the 
effective dynamics of QCD is described by a chiral Lagrangian whose target space is 
SU(N c xNf)L x SU(N c xNf)n, resulting from misalignment of the qq condensate in the 
color space. One then makes the additional assumption that the solitons of this effective 
theory are stable and may be quantized semiclassically. Remarkably enough, it turns out 
that such solitons have the same quantum numbers as the constituent quarks in the quark 
model. Thus the constituent quarks in this model are "skyrmions" in color space. 

This is a very attractive idea, but some crucial elements of the puzzle are still missing, 
most importantly a derivation of the effective dynamics from QCD. However, recent work on 
QCD in 1 + 1 dimensions ]TT[| (QCD 2 ) provides evidence in support of this picture. Exact non- 
abelian bosonization allows the Lagrangian of QCD 2 to be re-expressed in terms of a chiral 
field u(x) G U(N c xNf) and of the usual gauge field. This bosonized Lagrangian has topo- 
logically nontrivial static solutions that have the quantum numbers of baryons and mesons, 
constructed out of constituent quark and anti-quark solitons. The identification of these 
solitons with constituent, rather than current quarks, goes beyond the usual fermion-boson 
duality in two dimensions. It is based on the fact that their mass is generated dynamically 
and depends on the coupling constant. 

When N c solitons of this type are combined, a static, finite-energy, color singlet solution 
is formed, corresponding to a baryon. Similarly, static meson solutions are formed out 
of a soliton and an anti-soliton of different flavors. The stability of the mesons against 
annihilation is ensured by flavor conservation. These results can be viewed as a derivation 
of the constituent quark model in QCD 2 . Thus the idea of constituent quarks as solitons of 
a Lagrangian with colored chiral fields becomes exact in D = 1+1. 

Pending a more formal justification of the effective dynamics in D = 3+1, it is interesting 
to see whether the quark solitons or "qualitons" have phenomenologically acceptable static 
properties, beyond their quantum numbers. The model discussed in Ref. is the simplest 
one with the required properties. It is the analogue of the Skyrme model in color space. 

It is quite probable that the actual effective Lagrangian in four dimensions is consid- 
erably more complicated, but retains the same transformation properties under flavor and 
color. Past experience with skyrmion-like models shows that their qualitative predictions are 
remarkably independent of specific details of the Lagrangian jjTj| pL3| . The group theoretical 
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analysis of collective-coordinate quantization J?J is thus likely to be correct, regardless of the 
specific dynamics. So until the correct approximate dynamics is found, it is worthwhile to 
investigate in some detail the simplest model possessing the right symmetries. 



2 Model and Solution 



We start from the Lagrangian[|7| 
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plus a Wess-Zumino term, where is the gluon field strength, fi x is the chiral symmetry 
breaking scale, / is analogous to the pion decay constant f n , T a is a color SU(3) generator 
normalized by TrT a T b = \5 ab) R p = (eW^E) //j x , £ = exp(2TI a T a //), 0^ = 0^ + iA^T a is 
the color covariant derivative, a s = g^j^K is the effective color fine structure constant and 
the dots (• ■ ■) refer to other operators with four or more derivatives and higher powers of a s . 

Of the five parameters in ([]]) we have some experimental information only on / and a s , 
and even that is quite limited. The two-flavor model's axial flavor current suggests]?] that to 
leading order in a s , f ~ /„•/ v3 ~ 53.7 MeV, but higher derivative operators ignored in this 
treatment might renormalize / away from that value, just as in the usual Skyrme model 0. 
Here a s is the strength of the effective color interaction between the color-carrying condensate 
and the gluon field, in principle derivable by nonperturbatively matching the effective and 
fundamental theories 0. As no such calculation is available at present, in practice we can 
take as a rough estimate the value a s ~ 0.28 provided by the non-relativistic quark model. 
For the other parameters, dimensional analysis suggests [[[], 01 that /x x ~ 47r/ and c 4 ~ 1 ~ c a . 

Given a parameter set {a s , c 4 , c a , f, /i x }, we seek classical static solutions for both E and 
gauge fields, adopting a radially-symmetric skyrmion ansatz 
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where a c acts on a two-dimensional subspace of SU(3) C . To obtain the classical gluon 
configuration we note that it should have the same symmetry as S under combined spin 
and color isospin rotations. We work in A$ = gauge and eliminate the remaining gauge 
freedom, imposing the Gauss law constraint for unit winding number. We also note that 
energy minimization will remove the pure Yang-Mills terms, resulting in a configuration 
whose nonvanishing components are 
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Then the radial energy density of a qualiton (relative to the classical vacuum) is 
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and the classical contribution to the rest mass is 



drp(r) . 



(4) 



(5) 



The five terms represent the respective contributions of the color electric field, the color 
magnetic field, the kinetic term, the four-derivative term, and the symmetry-breaking term. 
Minimizing the classical mass yields the variational equations for F(r) and 7(r): 
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subject to the boundary conditions (Q) and 7(0) = = 7(00). The equations were solved 



with the COLSYS package fl4j . 

Given a particular solution {E,A}, minimizing the classical mass the symmetry 
generates a family of new solutions with the same energy: 



where f2(r, t) G SU(3) is a collective coordinate. Quantizing these zero modes in the gauged 
caseJ7| is more subtle than in the global-symmetry Skyrme model because the quantum states 
have to satisfy the Gauss law constraint. In the external background {£, A} corresponding 
to a given parameter set this constraint leads to differential equations for Q. We solve these 
equations numerically, obtaining the two moments of inertia I\ and I2, which govern the 
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contribution from quantization of the collective coordinates in the full semiclassical mass 
formula 

M = M d+ jU + i)(±-±-) + ±-(c 2 -±) , (9) 

where C 2 and j are respectively the SU(3) Casimir and the spin in the appropriate repre- 
sentation. The ground state solitons are color triplets with j = |, and mass 

M = M d + -|- + ±r = M d + M x + M 2 . (10) 

oli 4i 2 

In the analogous expression in the usual Skyrme model the quantum contribution, Mi + M 2 , 
is suppressed by 1/iV 2 relative to the classical mass. This justifies the semiclassical approach 
in the Skyrme model. In the present case, however, the large-iV c expansion does not help, 
because both the classical mass and the moments of inertia scale like ~ iV?: parametrically 
the classical and quantum terms are of the same order. For the semiclassical approach 
to make sense, the parameters of the model have to be chosen so as to make the quantum 
contribution numerically small, while at the same time giving phenomenologically acceptable 
results. 

Starting with the "natural" guess discussed earlier, {a s = 0.28, C4 = 1, c a — 1, / = 53.7 
MeV, // x = 4vr/}, we obtain M = M d +Mi+M 2 = (361+605+443) = 1409 MeV. This is 
doubly unacceptable: M is four times the phenomenological value of 350 MeV, and the ratio 
of quantum corrections to classical contribution is large, 77 = (M 1 +M 2 )/M c ; = 2.2. We can 
also compute the "half-height radius" ri defined by F(ri) = it/2, and obtain ri = 0.354 fm, 
which does look reasonable compared to the proton charge radius ~ 0.86 fm. T" his situation 
is rather generic: it is easy to find parameters yielding a smaller mass with a bigger radius or 
vice versa, but one cannot keep both of them small, while decreasing 77: moments of inertia 
scale like M d r 2 and so from flT0|) rj ~ l/(M c ^r) 2 . 

We can learn more by systematically investigating the parameter space. This task can 
be simplified by expressing all dimensionful parameters in terms of /. Then the XSB scale 
fM x appears only in the dimensionless combinations £4 = (47r///z x ) 2 C4 and c a = (yU x /47r/) 2 CQ,. 
Together with a s they form a three dimensional parameter space, which determines the 
physics up to an overall scale factor. A variational estimate then suggests that M has a 
local minimum, ~ 20/, near £4 ~ 3. In the c a direction the minimum occurs only at the 
boundary, c a = 0. Near the minimum the mass changes slowly as a function of 64, allowing 
partial minimization of either 7] ~ 2/(c4) 2 , or of rif ~ Q.ly/cl. 

Numerical results confirm these expectations. With a, = 0.28 and local mini- 

mum M = 18.3/ occurs at c 4 = 2.72, yielding 77 = 0.377. When c a is fine-tuned small, the 
minimum mass becomes M ~ 16/, falling to M « 13/ for a s = 1. Thus for any reasonable 
values of c a and a s , M is well above the colored pseudo-Goldstone mass mn = (|c Q 47ra s ) 1//2 /, 
and comparable to fi x ~ 47r/. Furthermore, a 350 MeV constituent quark requires / to be 
considerably below 53.7 MeV, implying a large value for ri ~ OSby/c^ (53.7 MeV//) fm. 
Clearly, an equally valid alternative is to take a larger /, so as to get a smaller ri, at the 
expense of increasing M above 350 MeV. For example, / = 36 MeV with c a small yields 
ri = 0.86 fm, M = 570 MeV, with rj unchanged. 
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At first it therefore seems that to obtain realistic properties of constituent quarks one 
must modify the model — which, as noted earlier, is only the simplest of a large family 
of Lagrangians possessing the right symmetries. The variational argument indicates that 
additional four-derivative terms are unlikely to make a qualitative difference, but it also 
indicates possibilities for quantitative improvements. 

We would like to point out, however, that there is another, somewhat unusual possibility: 
the radius of a constituent quark might actually be larger than the radius of the proton. We 
are used to composite objects being larger than the building blocks out of which they are 
constructed, but there are some exceptions to this rule when the interactions among the 
constituents are very strong. Clearly, when something like that happens, there must be a 
very substantial deformation of the original constituents. An example is provided by the 
fusion of two Deuterium nuclei|15| to form a nucleus of Helium-4. The r.m.s. charge radius 
of 2 H is 2.12 fm, with binding energy of 2.2 MeV, while the corresponding figures for 4 He 
are 1.67 fm and 28.3 MeV. The confining interaction between qualitons is stronger than the 
interactions that shape the qualiton. It is therefore possible that a similar phenomenon occurs 
when qualitons bind together to form hadrons. We do not know whether Nature realizes 
this scenario, or whether the effective Lagrangian then remains a useful guide, but if it does, 
we can find a broad region of the parameter space corresponding to a phenomenologically 
acceptable mass. 

To this effect, we first choose a somewhat small c a , to minimize the dimensionless mass. 
This will allow M = 350 MeV to correspond to a slightly larger /. Next we decrease C4 in 
the neighborhood of the minimum, so as to decrease the radius, while constraining rj < 0.33 
to protect the validity of the semi-classical approximation. With a s = 0.28, c a = 0.1, and 

ii y = iii f, we thus find parameters c 4 = 2.46,/ = 21.4 MeV, giving ri = 1.43 fm and 

2 

rj = 0.33. Fig. 1 plots the resulting shape functions. 

To test the sensitivity of this choice we note that taking instead c a = 0.3 decreases / 
only slightly to 20.6 MeV, and increases ri to 1.50 fm, with a range 1.44 fm < ri < 1.56 fm 
corresponding to 0.1 < a s < 1. Even taking c a — l,a s — 1 raises n only to 1.67 fm. Thus 
in the phenomenologically relevant regions of parameter space the results are rather stable 
and do not require fine-tuning. 



3 Summary and outlook 

We have studied the static properties of constituent quark solitons for the simplest case 
Nf—1, N c =3. Self-consistency of the semiclassical calculation indicates that in order to 
obtain phenomenologically acceptable constituent quark masses, such objects are likely to 
have rather large radii, before binding into hadrons. A more detailed investigation of this 
possibility requires extending the present work to multi-soliton configurations. Such an 
extension would also make it possible to compute modifications of constituent quark color 
magnetic moments in bound states, giving baryon hyperfine mass splittings and thus fixing 
the effective a s . 

It would also be interesting to extend the calculation to an Nf>2 model, since for Nf—1 
there is no chiral flavor symmetry and hence no pion. This would allow a computation of 
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the constituent quark axial- vector coupling g\. The framework can also be fruitfully applied 
in the study of the spin and strangeness content of the proton, or in a computation of the 
Isgur-Wise function [16] representing the light-quark degrees of freedom in the formalism of 
the heavy quark symmetry. 
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Figure Caption 

The profiles F(r), 7(7"), and the classical mass density p(r) in MeV/fm, with a s = 0.28, 
c 4 = 2.46, c a = 0.1, / = 21.4 MeV, and fM x = Anf. M = 350 MeV, 77 = M quantum /M cl = 0.33. 
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